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Based on density functional theory calculations and group theoretical analysis, we have studied 
NaLaMnWOe compound which has been recently synthesized [Phys. Rev. B 79, 224428 (2009)] 
and belongs to the AA' BB'Oe family of double perovskites. At low temperature, the structure 
has monoclinic P2\ symmetry, with layered ordering of the Na and La ions and rocksalt ordering 
of Mn and W ions. The Mn atoms show an antiferromagnetic (AFM) collinear spin ordering, and 
the compound has been reported as a potential multiferroic. By comparing the low symmetry 
structure with a parent phase of P4/nmm symmetry, two distortion modes are found dominant. 
They correspond to MnC>6 and WOg octahedron tilt modes, often found in many simple perovskites. 
While in the latter these common tilting instabilities yield non-polar phases, in NaLaMnWOe the 
additional presence of the A- A cation ordering is sufficient to make these rigid unit modes as a 
source of the ferroelectricity. Through a trilinear coupling with the two unstable tilting modes, a 
significant polar distortion is induced, although the system has no intrinsic polar instability. The 
calculated electric polarization resulting from this polar distortion is as large as ~ 16 pC/cm 2 . 
Despite its secondary character, this polarization is coupled with the dominant tilting modes and 
its switching is bound to produce the switching of one of two tilts, enhancing in this way a possible 
interaction with the magnetic ordering. The transformation of common non-polar purely steric 
instabilities into sources of ferroelectricity through a controlled modification of the parent structure, 
as done here by the cation ordering, is a phenomenon to be further explored. 



Perovskite oxides are one of the most interesting and 
studied classes of inorganic compounds. They show many 
interesting properties such as ferroelectricity, ferromag- 
netism, superconductivity, colossal magnetoresistance, 
etc. 0] The prototype structure is the simple ABO3, 
where the three-dimensional framework of corner-sharing 
BOq octahedra form 12- coordinate cavities occupied by 
the larger A cations. The possibility of A- or B- sub- 
stitution gives rise to a large compositional range with 
different properties and symmetries. A well studied case 
is the 50:50 substitution of B- site in the double per- 
ovskites A2BB'Os, where the B and B' cations order 
in a rock-salt fashion. 0, Q Analogously, ordering of the 
A-cations can also occur, showing a strong preference 
for a layered arrang ement.ii Finally, simultaneous A- 
and 5-site cation ordering can also occur in compounds 
of type AA'BB'O^ These compounds are interesting 
from a structural point of view since they exhibit rock- 
salt ordering of B-site cations and layered ordering of 
A-site cations. However, they are relatively unexplored, 



because this type of ordering is very rare. 

0, 1 Clearly, 

the presence of four different distinct cation sites greatly 
increases the possibility of functional design in this class 
of compounds. NaLaMnWOe that has been recently 
synthesized [6j belongs to the A and B ordered double 
perovskites. Mn and W ions are octahedrally coordi- 
nated to O ions and these octahedra are strongly tilted: 
the tilting is represented as a + a + c~ in terms of Glazer 
notation[9]. It exhibits a polar crystal structure with 
monoclinic P2i symmetry and an AFM spin configura- 
tion with Neel temperature Tn = 10 K [6-8]. P2\ is a 
polar group where electric polarization is allowed. There- 
fore, the magnetic ordering at low temperature, suggests 
a possible multiferroic behavior [sHsl , show ing both spon- 
taneous ferroelectric and magnetic order. 10h13| |. Typi- 
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cal multiferroics belong to the group of perovskite transi- 
tion metal oxides, and include rare-earth manganites and 
ferrites {e.g. TbMn0 3 , HoMn0 3 HoMn 2 5 , LuFe 2 4 ). 
Other examples are bismuth compounds BiFeOs and 
BiMnOa, non-oxides, such as BaNiF4, and spinel chalco- 
genides, e.g. ZnCr 2 Se4 14 -2^. In these compounds fer- 
roelectricity is either an independent phenomenon (like 
in BiFeOs), where ferroelectricity can be caused by an 
intrinsic polar instability, or its existence is an induced 
effect of the magnetic ordering (as in TbMnOa), where 



2 



ferroelectricity is a secondary effect, i.e. improper ferro- 
electricity. 

Very recently, new mechanisms for the coexistence 
and coupling of the two electric and magnetic order- 
ing have been investigated. An intriguing possibility is 
that a coupling of non-polar instabilities can give rise 
to a polarization, as first theoretically predicted in Ref. 
[29j |. Perez-Mato et al. studied SrBi2Nb20g compound 
and found that the ferroelectricity is due to the inter- 
play of several degrees of freedom, surprisingly not all 
of them associated to unstable or nearly- unstable modes. 
In particular, a coupling between polarization and two 
rotational modes has been described. [13] More recently, 
Bousquet et al. has demonstrated that ferroelectricity is 
produced by local rotational modes in a SrTiOa/PbTiOs 
superlattice.[3l| Following this work, Benedeck and Fen- 
nie proposed that magnetoelectric coupling, weak fer- 
romagnetism, and ferroelectricity can develop from the 
combination of two lattice rotations, neither of which pro- 
duces ferroelectric properties individually. [32[ In other 
words, separate lattice distortions can not only produce 
ferroelectricity, but also modify the magnetic order and 
favor magnetoelectricity. 33l- 35j 

In this paper, we will demonstrate that, by means of ab 
initio density functional theory (DFT) calculations and 
group theoretical analysis, a similar mechanism works in 
the case of NaLaMnWC>6 producing a relatively large ~ 
16 fiC/cm 2 polarization: two tilting instabilities are not 
enough to induce ferroelectric polarization by themself, 
but their coupling is sufficient to induce a polar instabil- 
ity and, therefore, ferroelectricity. To this end, we will 
show that the additional degree of freedom provided by 
the cation A-A! ordering has an important role. The 
mechanism explained in the article by Fennie et aZ. [3^| 
and termed "hybrid improper ferroelectricity" is basi- 
cally the same as the one discussed here, although for 
other type of systems. The calculated polarization is 
large here, due to the intrinsic softness of the polar mode 
in NaLaMnWOg- The condensation of a hard mode 
through a trilinear coupling with two unstable modes is 
in fact a quite common mechanism in perovskite-related 
compounds and other materials [3(| [36|, and can also 
be the origin of the simultaneous condensation of several 
distinct soft-modes, effectively acting as a single order 
parameter [33] 

Electronic structure calculations and structural op- 
timizations were performed by using the "Vienna Ab 
initio Simulation Package" (VASP) and Projector Aug- 
mented Wave (PAW) pseudopotentials.[38| The Perdew- 
Becke-Erzenhof (PBE) was employed for the exchange- 
correlation potential. [39( The DFT+U method was used 
to properly take into account correlation effects 40] using 
U —4 eV. We used experimental lattice constants and we 
relaxed the internal coordinates 0. The magnetic struc- 
ture is an AFM collinear spin configuration with prop- 
agation vector k = (0.5,0,0.5) (see Fig. [[]). A 2x4x1 
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FIG. 1: Magnetic unit cell of NaLaMnWOe with propagation 
vector k = (0.5, 0, 0.5) which doubles the nuclear unit cell 
along a and c directions. The orange arrows indicate the 
spins of Mn atoms. 



Monkhost-Pack fc-point grid in the Brillouin zone was 
used. The Berry phase approach 4l|, |42| was employed 
to calculate the electric polarization. 

A relaxed structure within the experimentally ob- 
served P2i symmetry was first calculated. The relaxed 
atomic positions are in good agreement with experimen- 
tal values, as shown in Table Q] Mn-0 and W-0 cor- 
ner sharing octahedra are strongly tilted and rotated due 
to the GdFeC>3-type distortion. It is well known that 
when the ionic radius of the A and A' cation is small 
(r Na = 1.22 A, r La = 1.24 A 43]), this distortion is en- 
hanced. In fact, in our case, the Goldschmidt tolerance 
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Na (2a) 


0.2514 


0.2289 


-0.0008 


0.2489 


0.2223 


0.0004 


La (2a) 


0.2556 


0.2755 


0.4984 


0.2574 


0.2845 


0.5003 


Mn (2a) 


0.7571 


0.2507 


0.2389 


0.7517 


0.2526 


0.2392 


W (2a) 


0.7581 


0.2494 


0.7652 


0.7547 


0.2451 


0.7664 


Oi (2a) 


0.5345 


0.5017 


0.6944 


0.5461 


0.5118 


0.6973 


2 (2a) 


0.5136 


0.5364 


0.3017 


0.5222 


0.5365 


0.3021 


Os (2a) 


-0.0205 


-0.0584 


0.2250 


-0.0362 


-0.0681 


0.2197 


4 (2a) 


-0.0669 


-0.0358 


0.7819 


-0.0692 


-0.0482 


0.7840 


5 (2a) 


0.8279 


0.2285 


0.5111 


0.8320 


0.2252 


0.5090 


Oe (2a) 


0.6684 


0.2565 


-0.0166 


0.6664 


0.2718 


-0.0127 



TABLE I: Experimental [6j and theoretical optimized atomic 
coordinates of NaLaMnWOe with P2\ crystal structure. The 
lattice constants are fixed to experimental values: a — 5.5717 
A, b = 5.5970 A, c = 8.0155 A, and (3 = 90.225°. 
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FIG. 2: (Color online) (a) Group-subgroup tree connecting 
the paraelectric PA/nmm phase to ferroelectric P2\ phase, 
(b), (c), and (d) are the pattern of atomic displacements for 
the Ig, IY, and T§ distortion mode, respectively. Yellow, 
green, purple, grey, and red spheres indicate Na, La, Mn, W, 
and O atoms, respectively. 



factor, i.e. ((r Na ,La) + r )/[v / 2((?'Mn,w) + ro)]=0.88, is 
smaller than 1. (Mn-O) and (W-O) bond lengths are 
2.126 and 1.953 A, respectively. The smaller (W-O) bond 
length is consistent with experimental results, where W 
ions are strongly displaced from the centers of the octa- 
hedra due to a second order Jahn- Teller effect. 0]. 

The tiltings present in NaLaMnWOe can be described 
using the Glazer notation Q, but it is preferable for 
our purposes to describe them within a comprehensive 
quantitative mode decomposition of the structure. 36 1 
The P2i structure of NaLaMnWC>6 can be considered as 
the distorted structure of a parent paraelectric PA/nmm 
phase, which is the highest symmetry structure com- 
patible with the mentioned order of the A and B sites. 
This idealized prototype structure can be found using the 
pseudosymmetry searching software PSEUDO.\a<\ Ba- 
sically, the structure is the ideal perovskite except for 
the ordering of the cations and there is no BOq and 
B'Oq octahedral tilting or rotation. Using the symme- 
try mode analysis program AMPLIMODES[Ab\, one can 
decompose the distortion relating the relaxed monoclinic 
structure with this prototypical structure into symmetry- 
adapted distortion modes. Three symmetry-breaking 
distortions can be distinguished with symmetries labelled 

as r+ r- 



and T 5 . Figure HJb, c and d) show the pat- 
tern of atomic displacements associated with these three 



different distortion modes. The 
modes can be seen to be essentially two tilting modes of 
the octahedron framework. They correspond to the so- 
called rigid unit modes (RUMs) [46j , that are frequently 
unstable in simple perovskites. Their usual labels when 
described with respect to a Pm3m perovskite are i?4+ 
and M3+ and would yield a distorted non-polar struc- 
ture of Prima (Pbnm) symmetry. Here, the structure 
of the modes is the same, but the cation ordering im- 
plies that the reference parent structure has a smaller 
tetragonal symmetry with a larger unit cell, so that they 
correspond to T modes for this larger tetragonal super- 
cell. The third distortion mode labelled and shown in 
Figure [U[d) is the polar mode frozen, which should be re- 
sponsible for the ionic spontaneous polarization present 
in this compound. The atomic displacements for this 
mode are mostly along the crystallographic b direction 
(polar axis). Na and La atoms displace in opposite di- 
rections by 0.055 and 0.108 A respectively, while Mn and 
W atoms remain essentially fixed. The maximum oxygen 
displacement along the polar axis is 0.136 A. 

The relative strength of each mode can be quantified 
by its amplitude, Q. [36| In our case, Q r + and Q r - are 

equal to 1.34 and 0.94 A, respectively. On the other hand, 
the polar distortion has much smaller amplitude of 
0.51 A. The relative size of symmetry-adapted mode am- 
plitudes offer a valuable clue to interpret the mechanism 
giving rise to the polarization. The fact that Q r + and 
Q r - are much larger than Q r -, suggests that the two 

1 5 

first distortion modes are the primary structural distor- 
tions with respect to the prototype phase, acting as order 
parameters, while the polar mode is a secondary induced 
distortion. Indeed, the and distortion modes 
lower the symmetry to P2i/m and P42i2, respectively 
(isotropy subgroups). Note that none of them is a polar 
subgroup, although the second one breaks the inversion 
symmetry. This means that the structure resulting from 
these two dominant distortions, which is very close to 
the real structure, is a non-polar structure. However, 
these two modes are sufficient to reduce the symmetry to 
P2\, which allows the presence of a secondary induced 
polar distortion. This can be understood considering the 
graph representing the group-subgroup relations connect- 
ing the high-symmetry PA/nmm to the observed P2i, as 
shown in FigfJJa). It can been seen that the intersection 
of the two mentioned isotropy subgroups is the observed 
P2i group, as it is their largest common subgroup. This 
means that the combination of these two tilting modes 
are able to induce a polar phase (although they don't pro- 
duce a polarization by themselves). It must be stressed 
that the original i?4+ and M3+ unstable modes in pure 
Pm3m perovskites preserve in all cases a non-polar sym- 
metry. It is only the fact that the cation ordering reduces 
the parent high symmetry to PA/nmm, which permits, in 
this compound, that two typical unstable non-polar per- 



ovskite tilting modes produce a polar structure, and, as a 
consequence, an improper ferroelectric. To confirm this 
scenario we have investigated the stability of each of the 
three distortion modes present in the relaxed structure 
of NaLaMnWC>6. We have calculated the energy varia- 
tion as a function of the amplitude of each mode, sepa- 
rately. The results are shown in Fig. [3] For the and 
IY modes, we have a characteristic double- well shape, 
clearly indicating that they are unstable modes. On the 
other hand, the total energy variation of T§ mode is a 
positive parabola, showing that the mode is stable. As 
shown in the inset of Fig|3tc), there are minima around 
Q =0.5 and -0.5 for the F^. In this sense, this mode 
is also unstable. However, the amplitude of the is 
0.5086 A and corresponds to Q = 1. Therefore, in or- 
der to achieve this amplitude, the distortion needs to 
combine with primary distortions (T§ and r7 modes). 

The presence of the T§ distortion with a quite signif- 
icant amplitude can be explained as an induced effect 
through a symmetry-allowed anharmonic trilinear cou- 
pling with the primary tilting non-polar distortion modes 
(r^ and rj") Q r -Q r+ Q r - . This simple coupling is suf- 
ficient to explain the presence of a non-zero amplitude 
Q r -rx(Q r +Q r -) in order to minimize the energy, despite 
the mode itself being essentially stable. More impor- 
tantly, this trilinear coupling will correlate the possible 
switching of polarization and tilting modes. For a fixed 
orientation of the polar axis, the number of domains is 
four and will be characterized by independent changes of 
sign of the two tilting modes, while the secondary polar 
mode correlates with them according to the proportion- 
ality law indicated above. Therefore, polarization switch- 
ing will imply necessarily a change of sign of one of the 
two rotational modes. This may be important to enhance 
the coupling of polarization with the magnetic ordering. 
Our results demonstrate that NaLaMnWOg can be con- 
sidered an improper ferroelectric. 

We now turn our attention to the estimate of the value 
of spontaneous polarization in NaLaMnWOg. For the 
calculation, the P2i/m phase resulting from the sole 
presence of the distortion is employed as a virtual 
paraelectric state and A is used as a scaling parameter 
of all ionic displacements, i.e. A = for the paraelectric 
state (P2i/m) and A = 1 for the sum of the two distor- 
tion components rr+r^ leading to the relaxed struc- 
ture (P2i). The negative P (A = —1) state is built 
by ionic displacements opposite to the positive P state, 
with respect to the A = structure. This corresponds 
to a correlated switch of the two distortion modes T± 
and T^, as discussed above. Figure [4] shows (a) the 
total energy difference and (b) the variation of polar- 
ization as a function of A. The total energy difference 
A£\-a =1 = E\ — E\—x shows convex and symmetrical 
shape between positive and negative P states, thus repro- 
ducing the bi-stable shape characteristic of ferroelectric 
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FIG. 3: (a), (b), and (c) are total energy variations as a 
function of the amplitudes of the Fg, Fj~, and FJ modes, 
respectively. The inset in (c) shows a zoom of Q values in the 
[-!,+!] range. 



materials. The energy barrier between the paraelectric 
and ferroelectric states A-&\ =0 -A=i — 890 meV/unit cell 
corresponding to 110 meV/formula unit (f.u.) is two 
times larger than in the YMn03 case, where A_B~60 
meV/ f.u. |47j] Spontaneous ferroelectric polarization AP 
calculated by the Berry phase method is ~ 16 /iC/cm 2 , 
as shown in FigHJb). This value is very close to the 
polarization value ~ 14 /^C/cm 2 obtained by the point 
charge model with nominal charges (i.e. +le for Na, 
+3e for La, +2e for Mn, +6e for W, and -2e for O). This 
means that electronic effects are are negligible and ionic 
displacements are the main contribution to ferroelectric 
polarization. The relatively large value of polarization 
is consistent with the softness of the polar mode shown 
in Fig. |3l We have also checked that, by reducing the 
on-site Coulomb interaction U to zero in our calculation, 
the polarization is not significantly affected. This sug- 
gests that Mn-3<i states do not play an important role in 
the ferroelectricity in NaLaMnWOg- 

Finally, we feel important to discuss two points, al- 
though they go beyond the purpose of the present study. 
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However, we believe that they may stimulate further 
studies in this interesting class of compounds. First, 
we have argued that the primary order parameters driv- 
ing the transition in this compound are the two rotation 
modes, and not the polar mode, which happens to also 
be unstable. For example, as in the Bousquet's work on 
PbTi03/SrTi03,[3lj one rotation mode completely sup- 
presses the other rotation mode, and it is required that 
the polar mode and a rotation mode be present in order 
for the second rotation mode to freeze in, i.e., the two 
rotation modes actually compete with each other and 
only one wins. In the work of Benedek and Fennie on 
Ca3Mn207.[33j the polar mode is hard and the two ro- 
tation modes are compatible with each other, i.e., true 
rotation driven ferroelectricity. In our case, we did not 
investigate whether the two rotation modes compete with 
each other or they are compatible: to properly address 
this question, a detailed study of the anharmonic cou- 
plings is needed. Second, a peculiar property arise by 
studying the magnetoelectric coupling in this compound. 
By symmetry considerations, not reported here, one can 
show that: i) a weak AFM magnetic moment M y should 
exist in this compound with the same propagation vec- 
tor of the M x and M z components; ii) the magnitude of 
the Mj, component should be proportional to the prod- 
uct of the amplitudes of the polar distortion and of the 
primary frozen magnetic mode in the plane xz. This 
means that one should be able to see that this y anti- 
ferromagnetic component switches by switching the sign 
of the polar mode and one could expect that it increases 
linearly with an increase of the amplitude of the polar 
mode. In this way, one can simulate the effect of an elec- 
tric held along y, and if the y-AF component responds 
linearly, this would correspond to a quite peculiar mag- 
netoelectric effect: namely an AF ordering responding 
linearly to an electric field. These considerations rein- 
force the idea of interesting physics in this novel class of 
compounds. 

In conclusion, we have studied a new multiferroic ma- 
terial, NaLaMnWOe , and shown that the compound is an 
improper ferroelectric. The estimated ferroelectric polar- 
ization is ~ 16 /iC/cm 2 , a relatively high value compared 
with ferroelectrics of this type. A negligible polar insta- 
bility does exist [see inset of Fig. [3] (c)], but the additional 
A cation layer ordering of this double perovskite is suffi- 
cient to make ferroelectrically active some tilting modes 
of the octahedra that in simple perovskites and in B- 
ordered double perovskites only give place to non-polar 
phases. Through a trilinear coupling with the two unsta- 
ble tilting modes, a significant polarization is induced. 
We hope that this study will stimulate further investi- 
gation of cation ordering as a tool to convert ubiquitous 
well-known steric non-polar instabilities into mechanisms 
for producing improper ferroelectrics, as well as new mul- 
tiferroics. 

The research leading to these results has received 
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FIG. 4: (a) Total energy difference AE\_\ =1 as a function of 
parameter A (in eV/unit cell), (b) Adiabatic switching path 
along monoclinic b direction between positive and negative 
ferroelectric states calculated by the Berry phase method (in 
/iC/cm 2 ). 
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